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For nonlinear sigma-models in the unitary symmetry class, the non-Unear target space can be 
parameterized with cubic polynomials. This choice of coordinates has been known previously as 
', the Dyson-Maleev parameterization for spin systems, and we show that it can be applied to a wide 

(<~*) ■ range of sigma-models. The practical use of this parameterization includes simplification of diagram- 

' matic calculations (in perturbative methods) and of algebraic manipulations (in non-perturbative 

, approaches). We illustrate the use and specific issues of the Dyson-Maleev parameterization with 

three examples: the Keldysh sigma-model for time-dependent random Hamiltonians, the supersym- 
metric sigma-model for random matrices, and the supersymmetric transfer-matrix technique for 
quasi-one-dimensional disordered wires. We demonstrate that nonlinear sigma-models of unitary- 
Uke symmetry classes C and B/D also admit the Dyson-Maleev parameterization. 



- 

X 



I. INTRODUCTION 
A. Nonlinear sigma-models 



^ ■ Extensive studies of disordered systems in last decades have identified the nonlinear sigma-model (NLSM) formalism 
as a universal tool for describing the low-energy properties of such systemsii^i^i^i^i^iiiSii^. The underlying idea of the 
NLSM approach is that the low-energy physics of disordered systems is determined by soft collective excitations 
(usually referred to as diffusons and cooperons). Therefore one can integrate out the high-energy degrees of freedom 
and end up with an effective theory which contains only soft collective modes. The resulting theory is a NLSM, which 
C ' is a field theory formulated in terms of a matrix field Q subject to a nonlinear constraint = 1. Depending on the 
, way of handling disorder averaging, the Q-matrix can act either in the supersymmetrio^ii^ii^i^, replica^ or Keldysh^iS 
(~| • spaces. Besides that, the Q-field can be a function of continuous space and/or time coordinates, and also have an 
O \ additional matrix structure due to an additional symmetry of the initial Hamiltonian. 
^ The NLSM action usually has a hydrodynamic form, containing only the lowest powers of Q, e.g., 

■ S[Q]^Tr{\A,Qf^BQ), (1) 

^ ' where Tr implies the full trace, involving integration over continuous coordinates (in the supersymmetric formalism, 
, the trace over the superspace is the supertrace). The action ([1]) determines the weight e~^^^^ in the functional integral 

■ which should be performed over the manifold 

(N ; Q = U-^AU, (2) 

' where 

in the "retarded-advanced" space, and U span an appropriate group of rotations. The matrices (or operators) A and 
-B in ([T]) must commute with A: [A, A] = = [B, A], so that Q = A is a saddle point of the action. 
}J] ' Depending on the symmetries of the disordered problem, the initial Hamiltonian may have additional symmetries, 
\ which translate into additional linear constraints on the matrix Q (realized in terms of extending its dimension, in 
combination with imposing constraints on the rotations 

The NLSM is a complicated field theory, and its exact solution is possible only in a few exceptional cases: zero- 
dimensional supersymmetric sigma-model for level statistics of random matrices^i^, and one-dimensional supersym- 
metric sigma-model for quasi-one-dimensional localization^i^iii. In a situation where the integration over the whole 
manifold of Q matrices cannot be performed exactly, the sigma-model can be treated with the standard perturbative 
approach, applicable in the weak-coupling limit (large conductance). Solution of the NLSM in the strong-coupling 
limit (small conductance) is a challenging and so far unresolved task. 

In the weak-coupling limit, a widely used perturbative expansion of the NLSM around the saddle point A is based 
on the parameterization of the curved Q-matrix manifold in terms of an unconstrained matrix W: 

Q^Kfm. ^=(_°5o)- (4) 
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where f{x) = 1 + x + + c^x'^ + (03 — l/8)a:^ + . . . can be an arbitrary function satisfying f{x)f{—x) = 1. 
A particular choice of the function f{x) is a matter of convenience (see discussion of various parameterizations in 
Ref. fl^ . and the calculated correlation functions do not depend on it. The integration contour in the space of b and 
b should be chosen consistently with the convergence of the integrals performed in the derivation of the NLSM. In 
the supersymmetric formalism, this amounts to the requirement of the compact fermionic and non-compact bosonic 
sectors^ii^. 

If the action ([T|) possesses several saddle points, an analogous perturbative expansion should be carried out around 
each of them. 

Substituting Q expressed in terms of W into the initial action ((T|), one arrives at a new theory 

J ... e-^I'^l DQ^ J ... e-^^[^] J/[iy] DW, (5) 

with the action S'/[W] ~ S[Af{W)] and, generally speaking, with a Jacobian J/[iy]. This new theory can be treated 
with the help of the standard perturbative approach by separating the action into the Gaussian part, Sj: [W], and 
the rest taken into account perturbatively. However, the initial action ([T]), which was a finite polynomial in terms of 
Q, typically becomes an infinite scries in terms of W. This makes higher-order perturbative calculations very involved 
as one has to take into account a rapidly growing number of diagrams made of various higher-order vertices oc W"'. 



B. Dyson— Maleev parameterization 

For nonlinear sigma-models in the unitary symmetry class, there exists a remarkable parameterization which does 
not belong to the infinite-series family ([¥]). The non-linear target space can be parameterized by cubic polynomials: 

fl-bb/2 b~bbb/A\ 
b -l + bb/2j- 

The asymmetry between the b and b matrices is the price one has to pay for finite number of terms in the parame- 
terization It can be shown by a direct calculation that the Jacobian of this transformation is unity, and thus the 
integration measure over b and b is flat: DQ = DbDb. In the new variables, the NLSM action becomes a finite-order 
polynomial, e.g., the action ([1]) becomes a sum of the bilinear and quartic in 6, b terms: 

s[Q\ = s^''\b;b] + s^''\b;b]. (7) 

The absence of higher-order interaction vertices greatly reduces the number of diagrams in the perturbative expansion, 
which considerably simplifies routine calculations. 

The parameterization ^ is closely related to the famous Dyson-Malee v^^i^^ (DM) parameterization for quantum 
spins. In that representation, the spin-S" operators are expressed by the boson creation and annihilation operators 
and d as 

S+ = {2S-a)a)a, S' = a\ S"" = S - a) a. (8) 

The Dyson-Maleev parameterization produces the correct spin commutation relations but violates the property 
(S'~)t = iS''^, and thus renders the spin Hamiltonian manifestly non-Hermitian. An alternative approach based 
on the Holstcin-Primakoff parameterizationi^ [which is in fact a counterpart of the parameterization ^ with 
f{W) = Vl + W^ + W] respects Hermiticity but generates an infinite series of interaction vertices. The Dyson-Maleev 
transformation has proven to be the most convenient tool for studying spin-wave interaction in ferromagneta^^^: it 
reproduces all the perturbative results obtained with the Holstein-Primakoff parameterization in a much faster and 
compact way. Owing to the analogy with the Dyson-Maleev representation of spin operators, the transformation Q 
will be referred as the Dyson-Maleev representation of the NLSM. 

In the original DM representation for spin operators, the Hilbert space of free bosons should be truncated in order 
for the operator S'" to have a bounded spectrum. This truncation of the bosonic Hilbert space is irrelevant for 
perturbative calculation but becomes essential in the non-pcrturbativc regime when the expectation value (d^d) is 
comparable to S. Clearly, some analogy of this Hilbert space truncation should be present also in the functional 
NLSM language. Indeed, in this paper we demonstrate that imposing certain conditions on the eigenvalues of the 
matrix bb makes the DM representation non-perturbatively equivalent to the original NLSM in terms of the Q field. 

At the perturbative level, the DM transformation ^ has been applied for the analysis of the replicatedi^ and 
Keldysbi^ sigma-models (cf. also the usage of the DM transformationi^ for two-dimensional classical ferromagnets). 
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At this level, when fluctuations are small, it suffices to use formally the parameterization ([6|) without taking care of the 
exact integration region over b and b. We illustrate the details of the perturbation theory in the DM parameterization 
with the examples of the Keldysh sigma- model for time-dependent random Hamiltonians in Sec. ITll 

At the non-perturbative level, the DM transformation ^ can be considered as an alternative parameterization of 
the complex manifold = 1. Therefore, changing variables from Q to b and b may be considered as a contour 
deformation on the NLSM manifold which leaves the functional integral invariant. We do not give a rigorous proof 
of this statement but demonstrate the non-pcrturbativc exactness of the DM transformation on the example of the 
supersymmetric NLSM for unitary random matrices in Sec. IIIII 

Another advantage of the DM representation is that it provides a compact algebraic way for deriving non- 
perturbative transfer-matrix equations for the one-dimensional NLSM. We demonstrate this with the diffusive su- 
persymmetric sigma-model for quasi-one-dimensional localization^ where the use of the DM transformation allows to 
derive the transfer-matrix Hamiltonian "without fermions" , using just the algebraic properties of the corresponding 
symmetric space, see Sec. ITVl 

Finally, in Sec. |V] we show that the DM transformation introduced for the unitary symmetry class can be ex- 
tended to other unitary-like symmetry classes C and B/D introduced by Altland and Zirnbauer— in the context of 
superconductivity. 

The main facts about the Dyson-Maleev transformation arc summarized in Sec. IVIl 



II. PERTURBATION THEORY: KELDYSH SIGMA-MODEL FOR TIME-DEPENDENT RANDOM 

HAMILTONIANS 



We illustrate the utility of the DM formalism in the perturbative regime with the problem of a quantum particle 
subject to a time-dependent random unitary Hamiltonian H(t). This problem has been studied in the framework of 
the Keldysh sigma-modeli^, and the use of the DM parameterization has been shown to greatly reduce the number 
of diagrams, in comparison with the usual infinite-series parameterization (j4|). 



A. The model 



In a semiclassical approximation, the dynamics of a particle subject to quantum evolution with a non-stationary 
Hamiltonian H{t) can be described as a diffusion process in the energy space. Diffusive spreading of the wave function 
is characterized by the diffusion coefficient D which determines the rate of the energy drift with time t: 

[E{t) - E{0)]^ = 2A^Dt, (9) 

with A being the mean level spacing. The dimensionless diffusion coefficient D depends on the rate of variation of 
the Hamiltonian measured by the dimensionless velocity of adiabatic energy levels: 



"-my- 

In the limit a 3> 1, the energy absorption is due to transitions in the continuous spectrum, and this regime can be 
described by the linear-response Kubo formula^ii^. In the opposite limit, a ^ 1, energy is absorbed during rare 
Landau-Zencr transitions between neighboring levels. In these limiting cases, the diffusion coefficient D{a) has been 
obtained by Wilkinson^l: 

f(/3/2)a, 

""^"^ = \c,«(/^+^)/^ a«l, ^''^ 

where /3 = 1,2, and 4 for the orthogonal, unitary, and symplectic ensembles, respectively, and are numerical 
coefficients. Quite surprisingly, for the unitary ensemble, C2 = 1 indicating that D(a) = a both in the limits of small 
and large a. 

Unfortunately, the methods used in Ref. [2l| to get the asymptotics (|1H) cannot be generalized to finite values of a. 
A general approach to calculating the full dependence D(a) based on the Keldysh sigma-model formalism has been 
derived in Ref. [23l . The field variable is the operator Q which is the integral kernel in time domain with values in 
2x2 matrices in the Keldysh space (analogous to the retarded-advanced space in the supersymmetric formalism). 
The Q-matrix is subject to the constraint = 1, where a convolution over time arguments is implied. As usual, the 



FIG. 1: An example of the 4-loop diagram for the diffuson self-energy (from Ref. [T^ . 



integration manifold is the orbit ([2]) of Ku' — Stt'<y3 (Pauli matrix in the Keldysh space) under unitary rotations U. 
In the dimensionless form (time measured in units of A^^), the Keldysh action for the linearly driven random unitary 
Hamiltonian takes the for m^^i^'^ 



S[Q] 



tr / dt{di - d2)Qtit2 



+ — tr 

ti=t2=t 4 



dtdt' {t-t'YQtt'Qt't 



(12) 



where the coupling parameter a is defined in Eq. pop . and tr denotes the trace over the two-dimensional Keldysh 
space. The parameter a plays the role of the dimensionless conductance which controls the strength of fluctuations 
around A: they are small for a 3> 1 and strong for a <C 1. 

The NLSM ([T^ contains full information about the function D{a) which can be expressed as^ 



1 

D{a) = — lim — — ^ 



(13) 



17=0 



where the diffuson I?,, {t) is defined through the correlation function of the off-diagonal components Qtft^ = ^'i^{<^^Qtit2 ) 
of the Q field: 



[DQ]e-'^'^^Qi:lQ[;l = ^Sit, 



to + is 



(14) 



B. Dyson-Maleev transformation 

In the DM parameterization ([6]), the Keldysh sigma-modcl action (fT2l) takes the form 
S[b, b] = l f fdh dt2 fei2 [{di + 82) + aih - ^2)'] ^21 - ^ [ dh dt2 dh dU {h - i2)(i3 - ^4) 612&23&34&41. (15) 



(16) 



The propagator of the quadratic part of the action (|15|) is given by 



where P^^-* (t) is the bare diffuson: 



(17) 



Substituting (jl7p into (fT3|) . one gets the large-a Kubo result D{a) = a. Treating the quartic term in the action p5|) 
as a perturbation, we can compute the perturbative series for the diffuson 'D^it) and, consequently, for D{a) in the 
limit of large a. From simple power counting^^, L-loop diagrams give a contribution to D{q) proportional to a^^^/'^. 
For the unitary ensemble, the number of loops must be even, which corresponds to expanding in powers of a^^/^: 
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In two loops, it has been found that ^2 = (the nuUification of this diagram is nontrivial, and we are not aware 
of any symmetry reasons for this fact)^'^. In the DM parameterization, there is just one two-loop diagram, while 
the standard infinite-series parameterization ([4]) generates two distinct diagrams, one made of a 6-order vertex, and 
another made of two 4-order vertices. While the DM parameterization appears to be the most compact way of getting 
the analytic expression in two loops, its use at this order of perturbation theory is not crucial. 

The coeflacient d4 is determined by four-loop diagrams. Its calculation performed in Ref. [13 was only possible with 
the use of the DM parameterization. Indeed, four-loop diagrams in an arbitrary infinite-series parameterization ^ 
would contain interaction vertices of the b and b fields up to the tenth order. Then even a classification of the four-loop 
diagrams becomes a sophisticated problem. On the other hand, in the DM formalism with the single quartic term in 
the action ([T5]), there are only 20 irreducible four-loop diagrams for the diffuson self-energy, see an example in Fig. [T] 
Some of them are related by symmetries, and the total number of nonequivalent diagrams is 9. They are given by 
seven-fold time integrals of the form 

where PeiTi, . . . , Tj) and 5*3 (Ti, . . . , Tj) arc homogeneous polynomials of degrees six and three, respectively. Numerical 
evaluation of these integrals shows that d4 is zero within the accuracy of the calculation and supplies the upper bound 
for its absolute value: |c?4| < 3 x 10"'*'. A further discussion of this result can be found in Ref. 

To conclude this Section, we mention that the same diagram classification applies to a wide class of NLSM with 
the action ([1]), whose Dyson-Malccv representation ([7]) contains only the Gaussian and quartic terms. In particular, 
quantum dynamics with an arbitrary time dependence of the Hamiltonian H{t) can be described by the Keldysh 
action ([12]) with the slightly modified last ter m^'^i^^ , thus also admitting the DM parameterization with the quartic 
interaction^^. Finally, various versions of diffusive sigma-models also belong to this class and have the same diagram 
classification in the DM representation. 




III. NON-PERTURBATIVE APPROACH: SUPERSYMMETRIC NLSM FOR UNITARY RANDOM 

MATRICES 

We start discussing the non-perturbativc aspects of the DM transformation with the simplest example of the 
supersymmetric NLSM for the spectral statistics of unitary random matrices. 

The theory is formulated in terms of the 4x4 supermatrix Q acting in the direct product of the superspace (BF) and 
the retarded-advanced space (RA). The pair correlation function of the density of states, = A'^{p{E)p{E + 

is given by the integral^i^ 

R(uj) = —Re[ (stikAQye-'^^^^DQ. (19) 
16 J 

Here A is the mean level spacing, A = diag(l, —1)ra: a-nd k = diag(l, —1)bf is the supersymmetry-breaking matrix. 
The supertrace is defined as the bosonic trace minus fermionic trace. The sigma-model action is 

S[Q] = -'-^ sir AQ. (20) 

With the non-compact BB sector and compact FF sector, the action (pO)) possesses two saddle points; the standard 
saddle point Q = A ("north pole of the fermionic sphere") and Q = Ak ("south pole of the fermionic sphere"), both 
belonging to the same connected component. The "south pole" is known to be responsible for oscillations in the pair 
correlation function^^. 

The DM transformation ^ may be considered as a special change of variables which parameterizes almost the 
whole complex manifold = 1. However the "south pole" Ak does not correspond to any finite b and b and can be 
achieved only as a limiting point for b = diag(0, 4:/p)bf, b = dia,g{0, p)bf at p 0. 

In Efetov's parameterization^il, the integration contour over Q is chosen in such a way that the eigenvalues A b and 
Ai? of the block Q^^ are real and bounded to the strip 

Ai3 > 1, -1<Af<1 (21) 

(the eigenvalues of the block Q^^ are the opposite: —As and —Xp). Now a transition from Efetov's parameterization 
of the Q-manifold to the DM parameterization can be realized as a deformation of the integration contour in a 
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multi-dimensional space. It will be convenient to impose an additional requirement that the eigenvalues of the block 
qRR (and the block Q^^ as well) do not change during the deformation. Since the eigenvalues of the product bb are 
given by 2(1 — As) and 2(1 — Af), this requirement constrains the eigenvalues of bb to be real and to belong to the 
corresponding strip. The constraint on the eigenvalues of bb is reminiscent of the truncation of the bosonic Hilbert 
space in the original DM transformation for spin operators (where it is formulated as the constraint on the eigenvalues 
of the number of bosons a^a). We may also note a similar restriction of the integration region in the theory of the 
DM parameterization for classical ferromagnetsiS. 

Throughout the deformation from Efetov's to DM parameterization, we may maintain the unchanged relation 
between the commuting components of b and b: 



fell bi2 

^21 622 



/611 


= 1 


(-bli 


bi2\ 


^621 


^■22) BP 


\ 621 


bh) 



(22) 



BF 

This choice of the real section of the complex manifold = 1 respects the correct domains for As and Xp- The 
independent elements are the two complex numbers bn and 622, and the four Grassmann numbers &12, &21, 612, &21 
(we do not use complex conjugation of Grassmann variables in our formalism). An arbitrary integral over Q can then 
be rewritten in the DM variables as 

Jj'[Q]DQ = J J'[Q{b,b)]e{\B-l)eil-\F)Oil + \F)DbDb, (23) 

where 9(x) is the step function, and Db Db is the standard flat measure: 

n; ^ Re &11 dim 611 dRc 622 dim &22 „ „ ,7 ,t 

DbDb= ^ d0i2 (3021 «0i2 ao2i. (24) 

The integral (|23p may be calculated straightforwardly in the variables 6y using the explicit expression for the eigen- 
values of the product bb: 

on \ ^ h h h , (&11612 + 6i2622)(62i&ii +622621) 

2(1 - As) = 611611 -I- 612621 + ? ; f ; = , (25) 

611611 - 622622 + 612621 - 621612 

2(1 - A,) = 622622 + 621612 + (b.A2+h2b22){b^i,+bM ^ (26) 

611611 - 622622 + 612621 - 621012 
These formulas allow a direct calculation of correlation functions in the DM parameterization. 



A. Technicalities of the integration over b and 6 



We illustrate the above approach by calculating the integral 

Z = J p{stikAQ)F{stTAQ)DQ = J p{stT kAQ)F{str AQ) e{XB - 1)6(1 - Xf)0{1 + XF)DbDb, (27) 

which for a suitable choice of the functions p and F gives the level- level correlation function . 

Note that in the flat parameterization used, str kAQ = 4 — (1/2) str A;(66 -1-66) = 4— 611611 — 622622 does not depend 
on Grassmann variables. A direct calculation (expanding in Grassmann variables) shows, quite surprisingly, that the 
same is true for two step functions in ([27]) : 9{Xb — 1) = (?(— 611611), 9{1 — A_f) = ^(622622)- Therefore, only two terms 
in the integrand ((27)) contain Grassmann variables: F(str AQ) = F(str66) and 9(1 + Xp)- Expanding these functions 
in Grassmann variables and evaluating the integral over them, after a straightforward but lengthy calculation we 
arrive at 

/OC />1 
dXB J dXp p"{2Xb + 2Xf)F{2Xb - 2Xf) + p{4:)F{Q), (28) 

where the second term is the anomaly contribution. In particular, it is responsible for the proper normalization 
Z = F{0) for p= 1. 

To find the level-level correlator p9|) . we should use ([28)) with p{x) = .x^/16 and F{x) ~ cxp{iTTujx/2A). Evaluating 
the integrals over A b and Af we arrive at the well-known result^l 



(29) 
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By reproducing the standard expression (|29p we demonstrate that restricting the integration region over b and b to 
the strip (|2ip makes the DM transformation non-perturbatively equivalent to the initial NLSM in the Q-representation. 

To be precise, our example p7|) involve only the diagonal blocks and Q^^. Therefore, our derivation is 

equivalent to that in the Hcrmitian "square-root-odd" parameterizatio n^^i^^ of the type @ with f{W) = 1-1- + 
W \J\ -I- which has unit Jacobian. This parameterization has the same diagonal blocks as the DM one and the 

same integration measure. 

In fact, a continuous interpolation is possible between the "square-root-odd" and DM parameterization, such that 
the diagonal blocks Q^^ and Q"^"^ remain unchanged, and the Jacobian equals one throughout the interpolations^,. 

The power of the DM parameterization will come in the full glory in the next Section when the off-diagonal blocks 
of the Q matrix will play the role. 



IV. NON-PERTURBATIVE APPROACH: SUPERSYMMETRIC NLSM FOR 
QUASI-ONE-DIMENSIONAL LOCALIZATION 

In this section, we apply the Dyson-Maleev parameterization to the transfer-matrix treatment of the Efetov su- 
persymmetric sigma-model for quasi-one-dimensional disordered systems. While most of the results of this section 
reproduce the known ones^*^, we find it instructive to rederive the transfer-matrix formulas with the DM method for 
two reasons: 

• the derivation of the transfer-matrix Hamiltonian simplifies in the DM parameterization and does not involve 
Grassmann variables, but uses algebraic relations for supermatrices instead; 

• the derivation of expressions for the correlation functions uses more transparently the symmetries of the sigma- 
model, and the intermediate states are more easily classified by the appropriate representations of the symmetry 
supergroup. 

The main goal of the calculation in this section is obtaining a closed expression for the two-point correlation function 
in the one-dimensional sigma-model^: 



[DQ] stT{AiQ{xi)) stT{A2Qix2)) e^^I^l , (30) 
where Ai and A2 are arbitrary supermatrices. The sigma-model action is 

S[Q] = -^sti J dx 



dx 



(31) 



where vi = vSw is the qlD density of states {v is the 3D density of states and Sw is the area of the wire), D is the 
diffusion coefficient, w is the energy difference between the retarded and advanced sectors. The supertrace str() is 
defined in the convention of Refs.|3|l5|: bosonic trace minus fermionic trace. We consider the unitary symmetry class, 
and Q is the supermatrix of dimension 2|2: it has the 2x2 structure in the retarded-advanced space. It is subject to 
the sigma-model constraint = 1 and is obtained by rotating the "north pole" ^ [independently at every position 

. . . 

First, we rederive the mapping of the path integral (j30p onto a finite-dimensional quantum mechanics (the transfer- 
matrix approach and then we express it in terms of dynamic correlation functions in this quantum mechanics. 



A. Transfer-matrix Hamiltonian. General formalism. 



It will be convenient to rescale the coordinate x in (|3ip by the localization length ^ = 2'kviD. The dimensionless 
length T ~ x/S^ will play the role of (imaginary) time in the quantum mechanics. In the rescaled variables, the action 



S = str / dr 



(32) 



where VL = -iuj^,'^ / {iD). 
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With the use of the DM parameterization ([6]), the action may be rewritten as: 

1 



S 



dr L . 



L 



■ str 



66- 



^6666 - 
4 



41766 



(33) 



This action can be converted into the transfer matrix (Hamiltonian) with the usual Lcgendre transformation. The 
momentum variables conjugate to b and b are 



n 



dL 
'db 



1- 



-bbb 



n 



dL 
di 



1, 



(34) 



(here 6, 6, H, and 11 are supermatrices, and their products involve the usual index convolutions), and the transfer 
matrix H is given by the Legendre transform 



H = L-stT 



(nh + nbj 



= A„i + Ao + C/, 



where 



A_i = 4 str (nn) , Ao = - str (UbUb) , U = -n str(66) , 



(35) 



(36) 



(Ao preserves the degree of 6 and 6, A_i lowers the degrees by one, and U raises the degrees by one). 

This Hamiltonian should now be understood as a quantum Hamiltonian with the canonical expressions for the 
momentum variables: 



H 



d_ 

db' 



H 



d_ 
db' 



(37) 



Since 6 and H do not commute, a care is needed for a proper ordering of the quartic term Ao of the Hamiltonian. 
However, in our example of the supersymmetric sigma-model, the commutators of 6 and H in the matrix product 
in pop generate terms proportional to str 1 = 0, therefore ordering is not important. The same situation occurs in 
the replica and Keldysh sigma-models, which are also defined so that tr 1 = 0. A more careful consideration may, 
however, be necessary if one wishes to extend this construction to other sigma-models (with non-normalized partition 
functions) . 

The above formal derivation is applicable to any sigma-model with the Dyson-Maleev parameterization. We further 
specify the explicit form of the Hamiltonian by 

• projecting the Hamiltonian onto a particular representation of the symmetry group of the sigma-model action 

dSID; 



• restricting 6 and b to belong to a particular finite-dimensional space. 



B. Hamiltonian in the invariant ("singlet") symmetry sector 



The sigma-model action (and hence the Hamiltonian H) is invariant with respect to the supergroup Hr x Ha of 
independent rotations in advanced and retarded sectors, as described in Appendix [XI In this subsection, we consider 
the restriction of the Hamiltonian ((35|) - ([36|) onto the sector of states invariant under the H]^ x Ha group (the "singlet" 
sector). 

In this paper we keep the symmetry analysis at a simple level and refer the interested reader to Ref. [sol . which 
contains a helpful introduction into the representation theory of the supergroup GL(1|1) relevant for our example. 

The "singlet" wave functions may generally be written as functions of cyclic traces x„ = str(66)". When projected 
onto such singlet functions, the Hamiltonian takes the form 



Ao^ A 



iS) 



E 

ki .k2- 



A- 



A 



4 E 

ki.k2 — l 



kik2 Xfci+fca-l 



_d d_ 

dXki dXk: 

d d 



+ (fcl + fe) Xk^Xk^ 



d 



dxk 



U 



dxki dxk_ 

[/(■^^ = -nxi 



+ (fci + A:2 + 1) Xkj^xk^ 



dx 



ki+k2 + l 



(38) 

(39) 
(40) 



Note that so far the formalism applies generally to arbitrary sigma-models admitting the DM parameterization. 
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C. Reducing the Hamiltonian to a finite-dimensional space 

Now we make the second step of specifying the particular type of the finite-dimensional symmetric (super)space. 
This will be done with the so-called "closing relations" . 

In our example, with the matrix Q having the dimension 2|2, the matrices b and h have the dimension 1|1, and 
so does their product hh. For this finite-dimensional space of matrices, among the infinite set of variables Xi, only 
the first two are independent. All the higher variables Xi may be expressed in terms of Xi and X2 with the help 
of "closing relations" which carry information about the properties of the particular symmetric superspace. When 
the Hamiltonian ([55)) . ([55)l - (pni) acts on a wave function of the form '^{xl,X2)^ only two higher-order variables are 
generated: x^ and x^. Thus it is sufhcient to use the two closing relations for the space of (1|1) supermatrices: 



2^3 = -p- [ixl + x\) , 
4x1 



2:4 = ^ 1 XIX2 + — 



(41) 



For future discussion of RA sector, it is helpful to note that those two relations are consequences of one matrix 
relation (the characteristic polynomial): 



X2 



(bby- -± bb + - 



Xl 



1 



0. 



(42) 



Using the closing relations (|4T|) . we reduce the operators ((38)) and ([39|l to the differential operators acting on the 
functions of two variables '^{xi,X2)- 



A, 



A, 



(s,i|i) 







X2 



d 
dxi 



A 



= 4 



Xi 



3^2 
Xi 

d 
dxi 



4X2 



_d d_ 

dxi dx2 

d d 



dxi dx2 V xi 



x\x2 



ixl 



X2 



x\ 



d 
dx2 

d 
dx2 



2x1 



d 

dx2 



(43) 
(44) 



Finally, it is convenient to change the variables to the eigenvalues \b and Af of the matrix Q^^ (or, equivalently, of the 
matrix —Q^^)- The matrix bb then has the eigenvalues 2(1 — As) and 2(1 — Af), which produces the changc-of-variable 
relations 



xi = 2(Af - As) , X2 = 4(Af - Ab)(2 - Af - As). 
With this change of variables, the Hamiltonian 

^(5,l|l) ^^(S,l|l)^^(S.l|l)^f;(5) 

is equivalent to that in Efetov's Eq. (11.38) in Ref. 0- 



(45) 



(46) 



D. Bilinear forms and self-conjugate singlet Hamiltonian 

The Hamiltonian derived above is not self-conjugate with respect to the flat integration in As and Af. For the 
convenience of interpreting it as a quantum-mechanical Hamiltonian, we will use an appropriate Jacobian to make it 
self-conjugate. 

For this step of the derivation, as well as for the future use, we will recall some well-known properties of the 
supersymmetric formalism. 

The "natural" bilinear form on wavcfunctions ^{Q) can be defined as 

(^-1(^)1*2(0))= j DQ^i{Q)^2{Q) = j DbDb^i^2, (47) 

where the integration is performed over the finite-dimensional superspace of Q with the normalization J DQ = 1. The 
second equality reflects the fact that the Jacobian of the DM parameterization equals one. The integration measure 
over DbDb is defined in dM]). 
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When restricted to the "singlet" states, this bilinear form is degenerate: it is given by the "anomaly"— 

(*l(AB,A;^)|«'2(AB,Af^)) =*i*2 . (48) 

As — Af — 1 

The anomaly gives 1 for the ground state and for the excited states. To define a nondegenerate bilinear form on 
the excited states, we can define the "renormalized" form 

(*i(Ab,Af)|*2(Ab,Af)), = J dXp d\B *i*2 J (49) 

with the Jacobian 

This bilinear form (-I-)* has the following properties: 

• with respect to this form, the Hamiltonian is self-conjugate. In fact, the Jacobian (j50ll can be obtained 
uniquely (up to an overall normalization) from this condition; 

• it is well-defined and non-degenerate on singlet excited states; 

• it is divergent on the ground state, and its divergence is proportional to (-I-). 
We will also use the "flat" bilinear form on singlet states: 

((*i(AB,Ai.)|«'2(AB,AF))) = J dXp dXs *i*2. (51) 

To produce a Hamiltonian self-conjugate with respect to this flat scalar product, it is sufficient to redefine the wave 
function by 

f{XBAF) = Vj^{XB,XF), VJ= , ^ , (52) 

Remarkably, in terms of the new wave functions /(AbjA^?), the Hamiltonian takes a very simple form, with the 
bosonic and fermionic degrees of freedom separated: 

Ho = VjH'-'^-^\'^^^HiXB)-HiXF) (53) 
V J 

where 

H{X)^-^{l-X')^ + 2nX (54) 

(the same expression in the bosonic and fermionic sectors). This expression is convenient to compare with Eq. (11.38) 
in Ref . i. 

The separation of variables implies that the excited eigenstates have the product form /b(Ab)/f(Af) with the 
energies 

E(b) _ E<-^), where E<-^^ and E'-''^ are the eigenvalues of the Hamiltonians H{Xb) and H{Xf), respectively 
(even though those two Hamiltonians have identical algebraic form, they are defined on different intervals Ab G [1; oo) 
and Xf G [—1; 1], and therefore have very different spectra). 

The separation of variables does not apply, however, to the ground state: it has the boundary condition ^'(Ab = 
Af = 1) = 1, and the corresponding wave function /(Af, Ab) given by (|52p is not normalizable. 



E. Separation of variables for the ground state 

It has been noted in Ref. [3l| that the ground state of the Hamiltonian ([53]) can also be obtained with a separation 
of variables, but with a different one (only applicable to the ground state). We would like therefore to make a short 
deviation from the main line of our paper (DM parameterization) to briefly review the construction of the ground 
state (in a slightly different way from the original derivation in Ref. Isil ) . 
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If one defines 



''--Tk-lk- '''' 



then a remarkable algebraic identity follows: 

Ho\/jD_ = D_y/JH* (56) 
where the variables separate again in the "dressed" Hamiltonian H^. 

H* ^H*{Xb)^H*{Xf), (57) 

with 

H*{X) = {1-X^)(^^^ +2nX. (58) 

With this algebraic trick, the zero-energy ground state of Hq is obtained from the zero-energy state of by applying 
the "dressing" operator \/jD^, 

fo = VJD^^^^HXB)^^P{XF), (59) 

where *^^^(Ab) and «'^^'(Af) are the bosonic and fermionic eigenstates of H'^{Xb) and H'^{Xp) at the same energy 
E. One can verify that the spectrum of H'^{Xb) is bounded by 2fJ from below and that that of H'^{Xf) is bounded 
by 2il from above, and therefore the only possibility to cancel the bosonic and fermionic energies in ([59|) is to 
take the bosonic and fermionic eigenstates at the same energy E = 2n. Those eigenstates arc given by modified 
Bessel functions, which results in the following explicit expression for the zero mode^ [with the proper normalization 
*o(As=Af = 1) = 1]: 

*° " 2ii7j^° " ^""^^^ ^B)PFh (pf)] , (60) 

where pB = 2^217(1 + As), Pf = 2^/20.(1 + Xf)- 

This alternative derivation of the ground-state wave function suggests the presence of a certain algebraic structure 
of the Hamiltonian related to the supersymmetry. Unfortunately, at present we are unable to properly identify this 
structure and leave this interesting question for future studies. 



F. Matrix-element reduction rules 



Our calculations in the DM parameterization may be simplified with the help of relations connecting the expectation 
values of non-singlet combinations of b and 6 to a singlet bilinear form in the coordinates Xb and Xf- These relations 
are somewhat similar to the Wigner-Eckart theorem, with the role of the symmetry group played by the supersymmetry 
Hji X Ha of the action. 

For the Q^^-Q^^ correlation functions (discussed in the following subsection), we will need the following relation: 

(*i|str(fcfl66) stT{kAbb)\^2) = 4 strfc^ strfc^ ((*i|*2)) (61) 

for any two Hfj x iJ^i-invariant states '^i{Xb, Xf) and 'I'2(As, Xf) and for any two supermatrices fc^ and kA- In this 
relation, the left-hand side involves the bilinear form ([47)) . and the right-hand side — the flat integration (jSTj) . 

This relation may be established by an explicit calculation (this is the only calculation in this section explicitly 
involving the Grassmann variables). In a somewhat different form, this relation has been used in Ref.Hllfor calculating 
local correlations (Q^^(O)Q'^^(O)). 

For the Q^^-Q^^ correlation functions (see discussion in subsection IIV H[) . we will need two other relations: 

str(fcflA6) str(fcAi?&)|*2) = str{kRAkAR) ((^-ilA/il^a)) (62) 

and 

str(fcflA6) stTikARbbb)\^2) = stiiknAkAR) ((*i|M2|*2» (63) 

for any two x iJ^-invariant states ^'i and and for any two supermatrices k^A and kAR- Here Mi and M2 are 
some functions of As and Xp. A direct calculation gives 

2 Ar + Ap — 2 , , 
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G. Diagonal two-point correlations 



We now turn to the calculation of the correlation functions (|30p . From symmetry considerations, one easily finds 
that the only nontrivial correlations arc those between RR-AA and between RA-AR blocks: 

(str[fci^g^'^^(0)]str[fcAQ^^(r)]) and (str[fcflAQ''-^(0)] strifc^iiQ-^^M]), (65) 

where k^, kA, k^A, and kAR are some arbitrary supermatrices. The RR-AA correlations will be further referred to 
as "diagonal" and RA-AR as "off-diagonal" . 

At T = 0, these correlations can be easily found from the explicit form of the ground state ^'o, together with the 
reduction rules discussed in the previous subsection. The results are known from Ref. Isi]. At t > 0, one needs to 
insert the appropriate evolution operator between the times and r. 

We will now specifically discuss the case of diagonal correlations. In this case, one of the two observables. Q^^, 
is invariant with respect to the Ha rotations, while the other observable, Q^^ is invariant with respect to the Hr 
rotations. Therefore one finds that all the intermediate excited states generated between and r should be singlets, 
and one can reduce the problem to calculating a certain expectation value of the singlet evolution operator exp(— iJo''')- 

A formal way to prove this statement is the following. First, we rewrite the RR-AA correlation function as 

{stT[kRQ'''^{0)]stT[kAQ^^{T)]) = (*o|str[fcAQ^^]exp(-i/r)str[fcflQ^«]|*o). (66) 
The state g^^^-o is a quadruplet of states of the form 

= IV'o + bhpi , (67) 



where tpn are some Hr x iJ^-invariant states. Under evolution with the Hamiltonian H given by ([35|) -([36 ll . this 
quadruplet evolves within a subspace spanned by the states (bb)"ipn. In our 1+1-dimensional example, the higher 
powers of bb can be expressed via the first two with the relation (|^^ . Therefore, the quadruplet exp(— i7T)(5^^^o 
can also be represented in the form (|67p . In principle, with the technique of the previous subsection, we can project 
the Hamiltonian onto such states and derive the corresponding differential 2x2 Hamiltonian acting on the functions 
{ipOjipi)- The projected Hamiltonian is obviously triangular: it preserves the subspace of Hr x iJ^-invariant states 
I'^o- In this subspace we can define a basis of eigenstates of H consisting of the ground state (with energy zero) 
and the excited states 

= stiibb)fkk' , fkk' = *f ^(As)*if^(Af ) (68) 
with the energies Ekk' = e\!^^ — E^fK We can complement this basis to the eigenstate basis in the space of states 



(j67| . Note that for any eigenstate quadruplet of the form (j67p . its supertrace should also be an eigenstate. Therefore, 
the eigenstates {ipo,tjji) have the form tpi = fkk', '4'a ~ 9kk' with fkk' given by (|68p and with some functions gkk' 
whose exact form is of no importance for our calculation. Now we can quite generally decompose 

Q^^*o = + ^^/^^'] + Mbb)fkk' + e^a (69) 

with some coefficients Ckk'-, dkk' , and e. By taking a convolution of this expression with (properly normalized) 5*0, we 
find e = 1. We further apply the evolution operator to this expression and take the convolution with Q^^'^q. Only 
the terms containing c^^/ and survive, and, using (|6ip . we find 



(str[fciiQ^^(0)]str[fcA(3'^'^(r)]) = -strfc^ str fc/j 



^Y.^lk'e-'''""^{fkk'\fkk'))+l 



kk' 



(70) 



Finally, observing that the coefficients Ckk' in (|69p may be obtained by decomposing in the basis fkk' of eigenstates 
of the Hamiltonian ([55]) - ([51)1 . we arrive at the final result of this calculation 

(str[fcHQ^«(0)]str[fc^g^^(r)]) = -strfc^ strfc^ [((*o|e-^«n*o)) + l] ■ (71) 



At T = 0, this expression agrees with the results of Ref. [sil 
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H. Off-diagonal correlations 

In a similar way, one can derive an expression for the off-diagonal two-point correlation function. The main difference 
is that now the intermediate excited states are of the form 

* = 6^0 + bbb<j)i (72) 

(with Hji X i/^-invariant (f>o and 0i), instead of (|67p . Both components c/jq and 0i are coupled to each other by 
the Hamiltonian, and thus we obtain an expression analogous to (|7ip . but with the effective Hamiltonian having an 
additional 2x2 matrix structure (in the (^^^(pi basis). 

From the same matrix closing relation (|42p . we can obtain the closing relations in the RA sector: 



X2_ 



bibbf = [±^]bbb+^l^b (73) 



\3 _ ^1 + "^■^2 , X2{Xi X2) J 



b{bby = '2 



4x{ Axf 



By using these closing relations, the Hamiltonian pS]) projected onto the states ([72]) may be represented as a 2x2 
matrix of differential operators acting on the vectors {(po, (f>i): 



^(flAi|i)^^(5,i|i) p \^+\ ^-..'^' hr^, (74) 




where iy^'^'i'i) is the Hamiltonian in the singlet sector multiplied by the unit 2x2 matrix. 

As in the singlet sector, we rewrite this Hamiltonian in the variables Xb and A;^ [using (|45p ] and rotate it to a 
self-conjugate form: 



2(Ab-1) 2{\F-l)\h^F 0\ _ /1 + A 



/As 


+ 1 


As 


- 1 



The resulting form of the Hamiltonian is 



(75) 



Hra = H,+ { _v,^^\^ (76) 



where Hq is given by ((53)) -((54 l) [the Hamiltonian in the singlet sector], and 



Vr = 2 ^""^ - \ , = 2^%HMZZi . (77) 

(As - XfY (As - XfY 

This Hamiltonian is self-conjugate with respect to the obvious extension of the flat bilinear form (j5ip to two-component 
vectors $ = (0o, Pi)'- 

/oo pi 
dXB J dXF Tr$f$2. (78) 

Proceeding along the same lines as in the singlet case, we write 

(str[fc^^Q^-^(0)] str[fc^flg^^(r)]) - (*o| str[fc^flg^^] cxpi^Hr) str[fc;^^Q«^] |*o) (79) 
and use the reduction rules ([62|l -(|63 p to arrive at the final result 

(str[fcfl,^g«^(0)] str[fcAfl,Q^^(T)]) = - str(fcfl,^fc^fl) ((*o ® exp(-i/HAT)|u ® *o)) (80) 

where 



« = «-M J,„)^-V(«) = ^( (81) 



-1/4; 4 \M2j Ab-AfV-VA|~T 

is an auxiliary two-component vector corresponding to the two-dimensional space (j72p . 

One can show^ that this formalism is equivalent to that in Eqs. (11.46) and (11.47) of Efetov's book^. 



14 



V. OTHER SYMMETRY CLASSES 

It does not appear possible to generalize the DM parameterization to the symplectic and orthogonal random-matrix 
symmetry classes: the corresponding Q matrices have additional constraints mixing the off-diagonal blocks Qba and 
Qabi and we cannot reconcile those constraints with the asymmetric DM parameterization in terms of b and b. 

However, a generalization is possible for the symmetry classes C and B/D. These are the "superconducting" sym- 
metry classes discussed by Altland and Zirnbauer— , with the unitary- type interlevel repulsion (/3 = 2). 

We specify our discussion below to the supersymmetric formulation of NLSM. 

For the symmetry class B/D, one may impose additional constraints on the matrices b and b: 

b^k = b, kb'^ = b, (82) 

where the transposition operation is understood in the supersymmetric sense (with changing the sign of one of the 
Grassmann components, see, e.g., Eq. (2.21) of Ref. 0), and k is the superparity operator, 

so that for any supermatix A, {A^)^ = kAk. The condition ((82|l used in the parameterization ([6]) implies 

Q = -70V\ ^=(-io)^^' 

which is the constraint on the Q matrix in the B/D symmetry classic. 
Similarly, the constraints 

b^k = -b, kb'^ = -b (85) 

produce the DM parameterization for the C symmetry class [the spaces of Q matrices in the classes C and B/D differ 
by the interchange of the fermionic and bosonic componentsiS, and thus ((85|) is obtained from ([82]) by replacing k by 

The constraints (|82p or (j85p reduce the number of independent components in b and 6 by a factor of two. 

Note that in the class B/D, the superspace of Q matrices contains two connected components, and the two subclasses 
B and D differ by the relative sign of the contributions from those components^. This sign choice coirrcsponds to the 
presence or absence of the zero mode in the original random- matrix theory. In the context of the DM parameterization, 
this implies that the parameterization ^ produces only one of the two connected components (containing the "north 
pole" A), and the second component (containing fcA) may be obtained by an appropriate non-local rotation (see, e.g., 
Ref. m for details). 

The only other symmetry class with /3 = 2 is the Chiral Unitary class (AIII in Ref. ITol ). So far we were not able to 
adapt the DM parameterization for this class (with the space of Q matrices of the A|A type). 

VI. SUMMARY 

In this paper we discussed the Dyson-Maleev parameterization of nonlinear sigma-models of the unitary symmetry 
class. Contrary to the commonly-used Hermitian infinite-series parameterizations Q, the non-Hcrmitian DM trans- 
formation ^ parameterizes the NLSM target space with finite-degree (cubic) polynomials. The Jacobian of the DM 
parameterization is unity. We find that the DM parameterization can be introduced only for the unitary symmetry 
class and its "relatives" C and B/D classes characterized by the same level-repulsion parameter /3 = 2. 

The main advantages of the DM parameterization include: 

• Simplification of the perturbative diagrammatic expansion. Due to the absence of higher-order interaction 
vertices in the DM representation, the number of diagrams is strongly reduced compared to the general infinite- 
scries parameterization ([5]) . Classification and identification of corresponding diagrams becomes a rather simple 
combinatorial task. 

• Possibility of obtaining non-perturbative results without resorting to Efetov-like parameterizations of the Q 
matrix. Instead, one can integrate directly over b and b with flat measure, with proper constraints on the 
eigenvalues of the product bb. 
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• An "algebraic" approach to the transfer-matrix treatment of one-dimensional diffusive NLSM. The derivation of 
the transfer matrix relies only on the symmetries of the action and on the algebraic structure of the superspace 
of Q matrices, without resorting to explicit manipulations with coordinates. Furthermore, the expressions for 
wave functions and correlations in terms of DM fields manifestly specify their symmetries, which simplifies 
manipulation, analysis, and interpretation of results. 

We expect that the Dyson-Maleev parameterization will be useful for perturbative treatment of various NLSM and, 
in particular, of the Finkelstein's replicated sigma-model for interacting systems^. The DM transformation might also 
be of importance for the study of two-dimensional localization in magnetic field since the topological term becomes 
quadratic in b and &, see Appendix[Bl Finally, the "algebraic" approach to the DM parameterization (using symmetries 
and algebraic properties instead of explicit coordinates) may provide a helpful tool for establishing relations between 
the three formulations of the NLSM: supcrsymmetric, replica, and Keldysh. 
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APPENDIX A: FINITE-DIMENSIONAL DM PARAMETERIZATIONS: SPACES AND SYMMETRIES 

In this Appendix, we describe the spaces and symmetries involved in the DM construction in its finite-dimensional 
version. We restrict ourselves to the supcrsymmetric formulation of the NLSM and to the case of the unitary symmetry 
class without additional constraints (see the discussion of constraints for symmetry classes C and B/D in Section fV)) . 
For the sake of generality, we take arbitrary dimensions of the bosonic and fermionic spaces: ms and m^?, respectively. 
This consideration will then be equally applicable to ordinary and super spaces. 

We start our construction with the two complex linear (super)spaces Lji and La (retarded and advanced sectors, 
respectively) of equal dimensions mB\'niF (wc take the dimensions equal for simplicity; the construction may be 
extended for Lr and La having different dimensions). Then we define h and h to be elements of the complex spaces 
of linear operators from La to and backwards: 

b: La-^ Lr, I-.Lr^ La- (A1) 

The products bb and bb are then linear maps Lr — > Lr and La La, respectively. 

The matrix Q in ([6]) is a linear operator acting in the space Lr © La- In the above example of Lr and La having 
equal dimensions m^lmi?, the matrix Q has the dimension 2mB|2m^. The condition = 1 defines a complex 
manifold in this linear space of matrices, and the DM parameterization allows us to parameterize its connected 
component containing the matrix ([3]). 

On this complex manifold of Q-matrices, wc consider the group generated by the invertible linear transformations 
in Lr and La {Hr = GL(Lr) and Ha = GL{La), respectively). The action of the group Hr x Ha on the Q-matrix 
^ may be written in terms of b and b as 

b^URbUX\ 

b^UAbU]^\ (A2) 

with Ur G Hr, Ua G Ha- 

Since the group Hr x Ha ~ GL(raB\mF) x GL{mB\mp) is the group of elements leaving A invariant, the complex 
(super)space of Q matrices (all possible rotations of A) may be described as 

GL{2mB\'2^raF) I {GL{mB\mF) x GL{mB\mF)\ (A3) 

This is the AIII|AIII symmetric space in the classification of Zirnbauer— , which corresponds to the unitary random- 
matrix theory (class A). 

Our illustrative supcrsymmetric examples in Sections IIIII and IIVI correspond to the case ijib = ^^^F = 1- In principle, 
the DM parameterization may also be used in non-supersymmetric models (with niB ^ ttif)'. for example, the simplest 
case niB = rriF = corresponds to the original DM parameterization ([5]) for spin S = 1/2. 

Another important group is the symmetry group of the NLSM action ([1]). In our examples in Sections IIIII and 
IIVI the symmetry group coincides with Hr x Ha (the stabilizer of A). Then the coordinate space for singlet wave 
functions (As, Xf) may be mathematically described as the double quotient 

[GL{mB\mF) x GL{mB\mF)] \GL{2mB\2mF)/ [GL{mB\mF) x GL{mB\mF)] ■ (A4) 
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APPENDIX B: TOPOLOGICAL TERM IN THE DM REPRESENTATION 

The sigma-model action describing the low-energy dynamics of the two-dimensional electron system subject to a 
strong perpendicular magnetic field is given by^ 

S=-^ J drtriVQf + ^ J dr tr e ^,QV ^QV ,Q , (Bl) 

where axx and axy are the mean-field longitudinal and Hall conductances respectively, and e^i, is the antisymmetric 
tensor. The last term in the action (|Bip is usually referred to as the topological term. 
In the Dyson-Maleev parameterization ([6|) , the action reads 

Note that the topological term becomes quadratic in b and 6, whereas all nonlinearity originates from the first term 
in the action (IBll). 
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